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Abstract 

In this paper, a new quantum mechanical method to extract energy 
from black holes with contracting horizons is proposed. The method 
is based on a gedanken experiment on quantum energy teleportation 
(QET), which has been recently proposed in quantum information 
theory. We consider this QET protocol for N massless fields in near- 
horizon regions of large-mass black holes with near-horizon geometry 
described by the Minkowski metric. For each field, a two- level spin is 
strongly coupled with the local quantum fluctuation outside the hori- 
zon during a short time period. After the measurement of N fields, N- 
bit information is obtained. During the measurement, positive-energy 
wave packets of the fields form and then fall into the black hole. The 
amount of excitation energy is independent of the measurement result. 
After absorption of the wave packets and increase of the black hole 
mass, a measurement-result-dependent local operation of the N fields 
is performed outside the horizon. Then, accompanying the extraction 
of positive energy from the quantum fluctuation by the operation, 
negative-energy wave packets of the fields form and then fall into the 
black hole, decreasing the black hole mass. This implies that a part of 
the absorbed positive energy emitted from the measurement devices is 
effectively retrieved from the black hole via the measurement results. 



1 Introduction 



Essentially, two methods are known by which energy can be extracted 
from black holes. One is the classical method that is based on the existence 
of ergospheres of rotating black holes. It includes the Penrose mechanism 
of test particles [I], the superradiance of an electromagnetic field, [2] and 
the Blandford-Znajek effect of accretion disks with magnetic field [3] for 
Kerr black holes. Moreover, it is possible to extract energy from Reissner- 
Nordstrom charged black holes because ergosphere-like structures of the black 
holes appear for charged particles [3]. Even though energy can be extracted 
from black holes using these classical methods, the horizon areas of the black 
holes do not decrease. If the quantum effect of matter is taken into account, 
energy extraction with a decrease in horizon area is possible. This is the 
second method of energy extraction, called the Hawking effect [5], which 
describes black hole evaporation. Black holes emit thermal radiation, ac- 
companied by generation of negative energy flux that falls into the black 
holes. The Hawking effect takes place spontaneously outside the horizon and 
cannot be stopped artificially. In this paper, a new quantum mechanical 
method of energy extraction from black holes is proposed. This method is 
based on a gedanken experiment on quantum energy teleportation (QET), 
which has been recently proposed in quantum information theory [B] — [10] - In 
this method, energy is extracted by a local quantum operation after a quan- 
tum measurement of quantum fluctuation of fields outside the horizons; thus, 
it is not a spontaneous phenomenon like the Hawking process. Besides, the 
time scale of the QET process is very short compared to the Hawking-process 
time scale of order of inverse of black- hole temperature. 

The protocols of QET consist of local operations and classical communi- 
cation. By measuring the local fluctuation induced by a zero-point oscillation 
in the ground state of a many-body quantum system and by announcing the 
measurement result to distant points, energy can be effectively teleported 
without breaking any physical laws including causality and local energy con- 
servation. The key point is that there exists quantum correlation between 
local fluctuations of different points in the ground state. Therefore, the mea- 
surement results of local fluctuation in a region include information about 
fluctuations in other regions. By selecting and performing a proper local 
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operation that is based on the announced information, the local zero-point 
oscillation of a place far from the measurement point can be squeezed, yield- 
ing negative energy density. During the local operation, surplus energy of 
the local fluctuation is emitted to the external systems and can be harnessed. 

We consider this QET protocol for N massless fields in the near-horizon 
regions of large-mass black holes with near-horizon geometry described by the 
Minkowski metric. The Hartle-Hawking state of the fields in the black hole 
spacetime, which is the thermal equilibrium state with the Hawking temper- 
ature [TTJ, is reduced to the Minkowski vacuum state in flat spacetime. For 
each field, a two-level spin is strongly coupled with the local quantum fluctu- 
ation outside the horizon during a short time period. The spin is measured 
to give 1-bit information about the fluctuation. After the measurement of 
N fields, iV-bit information, denoted by Dk (k — 1 ~ 2^), is obtained. The 
post-measurement quantum states outputting Dk are not the vacuum states 
but the states excited by the measurement devices. Positive-energy wave 
packets of the fields form and then fall into the black hole. The amount of 
excitation energy is independent of D^. After absorption of the wave packets 
and increase in the black hole mass, a Dfc-dependent local operation of the 
N fields is performed outside the horizon. Then, along with the extraction of 
positive energy from the fluctuation by the operation, negative-energy wave 
packets of the fields form and then fall into the black hole. The mass of the 
black hole decreases due to absorption of the negative flux. This means that a 
part of the absorbed positive energy emitted from the measurement devices is 
effectively retrieved from the black hole via the measurement results. Though 
energy extracted by QET does not exhibit thermal property like Hawking ra- 
diation, the energy extraction process involving measurement and the local 
operation, which generates a pair of excitations with positive and negative 
energy, is regarded as a controlled Hawking process. It is noteworthy that 
if a unitary operation corresponding to a different iV-bit sequence Dk'^k is 
performed, the amount of energy extracted by QET decreases. Without any 
information about the correct Dk, a unitary operation corresponding to a 
random iV-bit sequence does not extract energy from the fields, but gives en- 
ergy to the fields on average. This implies that a black hole, which swallows 
the wave packets generated by the measurement with output Dk-, remembers 
the information about Dk- Thus, black holes absorbing quantum states with 
different Dk are physically different even though their classical geometries 
are the same. 

The paper is organized as follows. In section 2, we briefly review the 
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horizon shift of large-mass black holes in classical theory, using the CGHS 
model. In section 3, we analyze a gedanken experiment on QET in a near- 
horizon region. The summary and discussion are provided in section 4. We 
adopt the natural unit c—h—1. 

2 Horizon Shift in Classical Gravity 



To briefly review the property of horizon shift in classical gravity, we adopt 
a two-dimensional soluble dilaton gravity model called the CGHS model [12] 
for simplicity. For higher-dimensional spherically symmetric black holes with 
large masses, we obtain essentially the same conclusion by treating falling 
matter with a spherically symmetric distribution. In the CGHS model, the 
action is given by 

where (p is a real scalar dilaton field, A is a positive constant, and f n denotes 
N real massless matter fields. The static black hole solutions are given by 

ds 2 = - (l - ye" 2A ^ dr 2 +(l- ye~ 2A ^ ' dr 2 , 

where M represents the black hole mass. The horizon radius is given by r# = 
2^hi(^). Using new coordinates defined by x ± = ± e2Ar ~ ^x e±Al ~i 

the metric is rewritten as 

o dx doc 

ds 1 = — — - — . 

1 — A 2 x + x~ 

The future horizon stays at x~ = in the new coordinates. Because the 
scalar curvature on the horizon is proportional to A 2 , a very small A yields 
the Minkowski spacetime as the near-horizon geometry. In the near-horizon 
region with \x^\ ^ 1/A, the metric is reduced to ds 2 = —dx + dx~. The 
matter field f n obeys the equation d + d-f n = in the region and can be 
decomposed into left- and right-moving solutions as follows: f n = fLn(x + ) + 
f'Rn{x~)- Solutions with falling matter into the black hole with mass M are 
also analytically obtained [12] as 
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o dx dec 
as = 

1 - A 2 x+x~ - ^ J x _l dy+ dz+T ++ (z+) ' 

where T ++ (x + ) is the incoming energy flux and is given by ^2 n=1 (d + fL n (x + )) 2 ■ 
Therefore, the metrics of general solutions in the x coordinates also become 
the flat metric ds 2 = —dx + dx~ when M ^> A and \x ± \ <C 1/A. Hence, 
the horizon shift due to absorption of matter cannot be observed in these 
coordinates of the near-horizon region. 

In order to observe the horizon shift, rescaled coordinates are available; 

these are given by X ± = \J~^x ± [12] ■ The metric is then rewritten as 
, 2 dX+dX- 

ds 2 = = . (1) 

f - A 2 X+X- - \ JZ dY+ /I dZ+T ++ { V /AZ+) 

From Eq. (pQ), the horizon shift is completed soon after the falling matter 
passes. If the support of T ++ (x + ) is given by (— oo,x+), the horizon stops 

soon after Xf = \J~^xt ■ The metric in the region with X + > is given 
by 

, , dX + dX~ 
ds 2 



M 1 
A 



A 2 X+ (X- - X H ) 



where M' is the enlarged black hole mass and the shifted future horizon 
stays at X~ = Xh = ~^\J^; fl^ T ++ (x + )dx + . Even for higher dimen- 
sional spherically symmetric black holes of Einstein gravity, the horizon also 
stabilizes soon after the falling matter passes if the matter takes the shape of 
spherical shells. This is because any spherically symmetric vacuum solution 
of the Einstein equation is static and equal to the Schwarzschild solution. 
After the horizon is shifted, the classical geometry of the black holes holds 
no information except the total energy. However, as explained in the next 
section, the black hole remembers more information quantum mechanically. 



3 Quantum Energy Teleportation 

The horizon of the initial black hole stays at x~ = 0. Let us consider N 
quantum fields /„ (n = 1 ~ N) in the near-horizon region with the x 
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coordinates. Because the fields are N copies of a massless scalar field, we 
can concentrate on one of the fields and disregard the index n for the time 
being. We consider a QET process with a time scale much shorter than 
inverse of the black hole mass which is the time scale of generation of the 
Hawking radiation. Thus, the Hawking process can be neglected in the later 
discussion. The equation of motion is written as <9 + <9_/ = 0, and the general 
solution is solved by the quantum left- and right-moving solutions as / = 
/+ + /_ (x~). The spacetime coordinates of the Minkowski metric are 
given by t — \ (x + + x~) and x = ~ (x + — x~). The canonical conjugate 
momentum operator of / (x) = f\ t =o is defined by fl(x) = d t f\ t=0 and it 



satisfies the standard commutation relation, 



f(x), U(x') 



iS (x — x') 



The left-moving wave / 4 
as 



can be expanded in terms of plane- wave modes 



o V4 



ITLU 



where (a^) is an annihilation (creation) operator of a left-moving particle 
and satisfies 



~ L ~L\ 
a ui a u>' 



5{u 



(2) 



The right-moving wave /_ (x ) can be also expanded in the same way using 
the plane-wave modes. The energy density operator is given by 
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E[X) 



U(xY : +: 



d x f(x] 



where :: denotes the normal order of creation- annihilation operators for the 
plain- wave modes. The Hamiltonian is given by H = e(x)dx. In 
the near-horizon region, the Hartle-Hawking state [11] is described by the 
Minkowski vacuum state |0), defined by H\0) = 0. It is also satisfied that 
(0\e(x) |0) = 0. Let us define the chiral momentum operators as 



U±(x)=IL(x)±d x f{x). 
Then, the energy density can be rewritten as 



e ix) 



x 



X 



(3) 
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Following the method given in [5J, we perform QET for the state |0) 
as follows. Figure 1 shows a spacetime diagram that expresses this QET 
experiment. Consider a probe system P of a two-level spin located in a 
small compact region [xi,^] satisfying x± > outside the horizon in order 
to detect fluctuations of /. In a similar way of particle-detector interaction 
of Unruh [TB], we introduce a measurement Hamiltonian between / and the 
spin probe such that 

H m {t) = g{t)G+®a y , 
where g(t) is a time-dependent real coupling constant, G + is given by 

r-oo 

G + = - + / A(x)fl + (x) dx, (4) 

^ J — CO 

X(x) is a real function with support [xi,x 2 ]> and a y is the ^-component of 
the Pauli matrices of the probe spin. We assume that the initial state of the 
probe is the up state | + |) of the z-component a z . The time dependence 
of g(t) is assumed to be generated by a device with an external energy sup- 
ply from spatial infinity, respecting local energy conservation. In the later 
analysis, we choose a sudden switching form such that g(t) = 5(t — 0). After 
the interaction is switched off, we measure the z-component a z for the probe 
spin. If the up or down state, +|) or — |), of a z is observed, we assign the 
bit number b = or 1, respectively, to the measurement result. The mea- 
surement is completed at t = +0. The time evolution of this measurement 
process with output b can be described by the measurement operators M b , 
which are often used in quantum information theory [13] and which satisfy 



M hP M\ 



Tr 

p 



VIp 



ft 



where p is an arbitrary density operator of the field, the time evolution 



operator V = Texp 



-i f_ H m {t')dt' 



is computed as exp 



-iG A 



0Y 



for 



t > 0, and the trace Trp is taken to the probe system. The measurement 
operators M b are evaluated as 



M h 



-I) 



.1 



exp 



— iG+ <g> a v 



Hence, we obtain the explicit expression of M b such that 
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M = cosG + , (5) 
Mi = sinG+. (6) 

For the vacuum state |0), the probability of obtaining b in the measurement 
(0|MjM 6 |0) is independent of b and is given by 1/2. The explicit calculation 
of (0\MlM b \0) is as follows. 



(o|m;m 6 |o) = (0| 



-If 



1 

2 ' 2 



cos 2G 



10} 



(0|sin 2 / X(x)U + (x)dx |0) 



1 
2 
1 

2" 



In the above calculation, sin A(x)fl + (x) dxj is calculated as the sum 

of the products of odd-number field operators and the correlation functions 
with odd- number field operators vanish for the vacuum state |0) of this free 
field. The post-measurement states of / for the result b are calculated as 



|V*> = V2M b \0) = ±= (e-f £ |A) 



A) 



where | ± A) are left-moving coherent states defined by 



± A) = exp 



x 2 



±z / A(x)II + (x) dx 



|0). 



(7) 



(8) 



The two states |"0o) an d IV'i) are non-orthogonal to each other with (V'olV'i) = 
(A| — A) 7^ because the measurement is not projective [13]. For convenience, 
let us introduce the Fourier transformation of X(x) as follows. 



A» 



X(x)e iu ' x dx, 



(9) 



where the coefficient satisfies A*(— u) = X(u). Note that the left-moving 
momentum operator n + can be expanded as 
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n+ 



dcu 



Using Eqs. ©, ©, an d (jit)]) , we obtain the following relation: 



exp 



■l l 2 



A(x)Il + (x) dx 



a£ exp 



X' 2 



X(x)fl + (x) dx 



Therefore, the coherent state in Eq. (jBJ) satisfies 



a£|±A> = TA/-A(w)|±A>, 

7T 



(10) 

7T 

(11) 

(12) 



which is derived from Eq. (TTTj) and a^|0) = 0. Using a^, the state | ± A) 
can be expressed as 



du—\X(u)\' 



exp 


H 







I ± A) = exp 

From this expression, we get the following inner products: 



|0>. 



(A|-A) = (-A|A) = (0|2A) =exp 



du-\\{uj)f 

7T 



(13) 



The expectational value of the Heisenberg operator of energy density 
e (x, t) for the post-measurement state in Eq. (JTj) is given by 



(i/> b \e (x, t) = - [(X\e (x, t) | A) + (-A|e (x, t) \ - A)] 
- (-l) 6 Im(A|e(x,t)| - A). 



(14) 



From Eq. f|T2|) . it is easily checked that lm.{\\e(x, t)\ — A) vanishes. In fact, 
we can show 



Im(A|e(x,t) | - A) 



Im 



duo 



UJ 



71 



1 2 



\(tu)e- iuJX + A(w)V 



0. 
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Besides, Eqs. (EJ), ( 1T2"j) . and (TTO]) yield the following relations: 



(A|e(a;,t)|A) = (— A|e (x, t) | — A) = (d+A(x+)) 2 . 

Substituting the above equations into Eq. (tHj), we can conclude that the 
measurement excites wave packets propagating to the horizon with an energy 
density that does not depend on b such that 

(i/> b \e(x,t)\ik)= (<9+A(x + )) 2 . (15) 

This result independent of the measurement result b is striking as compared 
to a particle detector model in the Rindler spacetime discussed by Unruh and 
Wald. In the reference |T7], a particle detector coupled with a thermal bath 
of a quantum field is analyzed which excites the field with positive energy 
when the detector makes a transition and is observed in a excited state. If 
the state of the detector is observed in the initial ground state, we have 
different excitation of the field. Especially, in a causally disconnected region, 
negative energy is induced by the detector. Thus amount of energy induced 
by the detector in the model has explicit dependence of the measurement 
results. However, in our POVM measurement model, the probe spin as a 
detector excites the quantum field independent of the measurement result as 
seen in Eq. f lT5|) . Moreover, soon after the measurement, the energy density 
for each post-measurement state vanishes except around the measurement 
point. The amount of total average excitation energy is evaluated as 

/oo 
(d + X(x)) 2 dx. (16) 
■oo 

The energy flux falls into the horizon at x~ = 0. Let us assume that at 
t = T, the left-moving excitation has already been eliminated from [xi,^] 
and satisfies X(x + T) = for x G [a^,^]- The average state at time T is 
expressed as 

PM = ^e-^M,|0)(0|M fe V^. 

b 

It is worth noting that the state pm is a strictly localized state defined by 
Knight [14J, because pM is locally the same as |0}(0| at t — T and satisfies 
Tr [pM£ (%)) — for x G [a^i,^]. Moreover, for each state \ipb), 2n-point 
functions of the left-moving operators 11+ with integer n are the same as 
those of |0) outside the horizon after time T. Let us consider a 2n-point 
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function (^|f[ + (xi + T) ■ ■ •n + (x 2 „ + T)\ip b ) and assume that all Xk satisfy 
Xk > and that they do not coincide with each other. Then, using Eq. (J7J), 
the function is given by 

(^ b \U + (x 1 +T)--'U + (x 2n + T)\i/j b ) 

= 1 T (A|fl + (a; 1 +T)---fl + (a; 2n + T)|A) 
2 [ +(-\\IL+(x 1 +T)---IL + (x 2n + T)\-\) 

- (-l) b Im(A|fl + (a; 1 + T) • ■ ■ fl + (x 2n +T)\— A). (17) 

Here, we have used the relation 

(X\h + (x 1 +T)---fl + (x 2n + T)\-X)* 
= (-A|n + (x 2n + T)---n + (x 1 +T)|A) 
= (-A|fl + (x 1 +T)---fl + (x 2 „ + T)|A), 

which is ensured by the commutation relation 

TL+(x), IL + (y)] =2i8'(x-y), (18) 

and the fact that all Xk + T do not coincide with each other. Because / 
is a free field, we can adopt the Wick's theorem for the composite operator 
U+(xi + T) ■ ■ ■ Tl + (x 2n + T) as 

= ^:fl + (a; Sl +T)---fl + (a; S2fc +T): 

n(0|n+(^i +T)fl + (x W2 +T)|0), (19) 



x 



where the contraction sum runs over only those contractions between opera- 
tors in different normal-ordered products. Then, using Eq. ffl9|) and the fact 
that \(xk + T) = 0, we get 

(±A|n+(xi + t) ■ • ■ n + (x 2n + t)| ± A) = (o\n+(xi + t) ■ ■ ■ ii + {x 2n + r)|o). 

. (20) 

This is because the normal ordered contributions (±A| : n + (x Sl +T) • • • Il + (x S2k + 
T) : | ± A) vanish: 

(±A| : fl + (x si + T) ■ ■ ■ tl + (x S2k + T) : | ± A) 
= 2 2k d x X(x Sl + T) ■ ■ ■ d x X(x S2k + T) = 0. 
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Thus, (±A|li + (xi +T) • • • n + (x 2n + T)| ±A) is reduced into a sum of products 
of n two-point functions for the ground state |0) and equal to (0|n + (xi + 
T)---n + (x 2n + T)|0>: 



(±A|n+(xi + T) • • • fl + (x 2n + T)| ± A) 

= ^(o|n + (x Wl + T)fi + (x W2 + t)|o) • ■ ■ (o|n 
= (o|rV(xi + T)---n + (z 2n + T)|o). 

By use of Eq. ffl9l) . it can be also proven that 

Im(A|n + (xi + T)---fl + (x 2n + T)| - A) = 

This is because the following relation holds: 

(X\:fl + (x Sl +T)---fl + (x S2k +T):\-X) 
= (-l) k Q( Xsi+ T)...Q(x S2k +T)(X\ 

where Q(x + T) is a real function defined by 



i + T)IL 



+ T)|0) 



0. 



(21) 



A), 



Q{x + T) 



du- 
ll L 



X(cu)e' 



-iu)(x+T) 



+ X(u 



,*Juj(x+T) 



and (A | — A) is also real as seen in Eq. ffl3|) . Moreover, taking into account 
the fact that all Xk do not coincide with each other, it can be observed that 
the two-point functions (0|ll + (x u , 1 + T)Il + (x W2 + T)\0) are real. Thus, by 
using the above results and Eq. (fl9|) . we obtain Eq. (|2T|) . By substitution 
of Eqs. (120]) and (I2T]) into Eq. (fl7|) . it is proven that the 2n-point function 
in the outside region after time T is the same as that of the vacuum state: 



(^ b \fl + {x 1 + T)---fl + {x 2n + T)\^ b ) 
= (0\fl + (x l + T)---fl + (x 2n + T)\0). 

It is worth noting that many-point functions of the energy density operator 
e(x, T) can be obtained from these 2n-point functions with the point-splitting 
regularization scheme [15J. Besides the left-moving mode, the many-point 
functions of the right-moving mode for l^) are the same as those for |0). 
Therefore, after the wave packets are absorbed by the black holes, the fluc- 
tuations in energy density for the state l^) are the same as those of the 
vacuum |0): 
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W> b \e(x u T) ■ ■ ■ e(x n , T) = (0\e(x u T) ■ ■ ■ e(x n , T) |0>. 

Because the massless field can affect the spacetimes of black holes only via 
energy density, the difference in the measurement result b yields no difference 
in black-hole spacetimes. The difference in b appears only in odd-number 
many-point functions of 11+ . For example, the one-point function of Il + for 
\ipb) is dependent on b and is given by 

f°° hi r i 
(4> b \tl + (x + T)\4> b ) = (-if / dw- \{u)e- iulix+T) + \{ooYe w{x+T) . (22) 

Jo 71 1 J 

We recall here that we have N copies of the field /. Thus, the measure- 
ment output bit (b n = or 1) is obtained for each field /„ and the total data 
of the measurements are recorded as Dk = (6162 • • • &jv) with k — 1 ~ 2^. 
If we take a large N, quantum fluctuation in the total energy density of 
the N fields is severely suppressed due to the law of large numbers, and 
the back reaction of the quantum fields to the black hole can be treated 
by a classical gravity equation of the metric. In the equation, the average 
energy flux as a source term is given by N(d + X(x + )) 2 . The falling flux 
yields a horizon shift of SX^ independent of b in the X ± coordinates as 

Figure 2 shows the horizon shift. The initial horizon (X - = 0) is denoted 
by Hi in the figure. The new horizon after the shift is denoted by Hp. It 
should be noted that a spacetime region becomes static soon after the flux 
passes through the region. If we do nothing more, Hp is the final event 
horizon. However, as seen below, if the unitary operation dependent on D k 
is performed, negative energy flux is created and the horizon shifts again. 
The receding horizon is denoted by Hqet- The argument does not change at 
all if the falling energy flux interacts with collapsing matter within the black 
hole. In Figure 2, right-moving matter represented by dashed lines collides 
with the falling positive-energy flux inside the horizon. 

At t — T, we perform a unitary operation on the quantum field /; the 
unitary operation is dependent on b and is given by 



Ub = exp 



/oo 
p(x)fl + (x) dx 
-00 



(23) 



where 6 is a real parameter fixed below and p(x) is a real function with the 
same support of X(x), [x±, x 2 ]- We do not require any classical communication 
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phase for this black-hole QET because the operation in Eq. (1231) is performed 
in the same region [x\, X2] as that where the measurement is executed. After 
the operation, the average state of the field / is given by 

p F = U b e- lTti M b \0}(0\Mle^Ul 
6=0,1 

Let us introduce an energy operator localized around the region [xi, X2] such 
that H B = JZ w (x) e (x) dx. Here, w (x) is a real window function with 
w(x) = 1 for x G [xi,x 2 ] and it rapidly decreases outside the region. The 
average amount of energy around the region is evaluated as 

E B = Tr [p F H B 

= J2( Q \Ml^ T6 UlH B U b e- iT6 M b \0). 
6=0,1 

In order to simplify the expression of Eg, a useful relation is available: 
UlH B U h = H B - 0(-l) b / d x p(x)fl+(x)dx + 9 2 {d x p(x)f dx. 

J XI J XI 

This is derived from Eqs. ffl8|) and fl23l . Thus, E B is rewritten as 



E B = J2mle^H B e-^M b \0) 

b 

-6^2(-l) b {Q\Ml / d xP (x)fl+(x + T)dxM b \0) 
b J X1 

+ 2 [ {d x p{x)fdx. 



(24) 



Because we assume that T > \x% — X\\, the following two relations can be 
directly proven: 

0, (25) 



e tTH H B e~ tTH , M b 



d x p(x)fl + (x + T)dx, M b 



0. 



(26) 
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From Eqs. (jSJ) and we are able to check the following relations: 

Y,MlM h = l, (27) 



^(-1) 6 M^M 6 = cos 



(28) 



Substituting Eqs. f l25|) -f l28|) into Eq. f )24|) . we obtain the following expression 
of E B : 



E B = {Q\j TH H B e-* 1 "\0) 



TH I 



■f2 



0/ dar9 x p(z) (0|n + (ar + T) cos 2^+1 |0) 



+ e 2 / (a x p(x)) 2 ^. 



Note that 



(0\e lTH H B e- iTH \0) = (0|^ B |0) 



iw(ar) (0|e(a;) \0)dx = 0. 



Therefore, the average energy E B = Tr 
puted as 

E B = -07? + fl 2 £ 
where £ = J^ 2 (^^(x)) 2 cfo and 



77 = / cted x p(x)(0|fl + (x + T) cos ( 2G+ ) |0). 



after the operation is com- 



(29) 



It is possible to simplify this expression of 77. Using Eqs. (1TU1) and (|T2|) and 
the relation given by 



cos(2G + )|0) = -(| + 2A)-|-2A» 



we obtain the following relation: 



(0|n + (x + T)cos [2G+) |0) = — (0|2A) / dujuj\(uj)e 

71 Jo 



-iui(x+T) 
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Substituting Eq. §§§ into the above relation, we obtain 




2 r°° 1 

= — 2A / dyX(y)- — 

* Jo {x-y + T) 



(30) 



By substituting Eq. (130]) into Eq. (1291 . we get the final expression of r\ as 



By fixing the parameter 9 such that 



so as to minimize Eb, it is proven that the average energy around [xi,:^] of 



By virtue of local energy conservation, this result implies that the unitary 
operation extracts positive energy +|-Eb| from the field to the external sys- 
tems. Taking into account the existence of N fields, the total amount of 
extracted energy is N\Eb\- The measurement step for = (101 •■•0) is 
shown in Figure 3. In Figure 4, the extraction step is depicted. This process 
is analogous to the standard Hawking process of black holes, that is, pair 
creation of particles with positive and negative energy outside the horizon 
[5]. Just like in the Hawking process, the created negative energy — N\Eb\ 
of the fields falls into the horizon and decreases the black hole mass, accom- 
panying the second shift of the horizon (Hp — > Hqet), as shown in Figure 2. 
For each field, the sum of the positive energy generated by the measurement 
and the negative energy generated by the operation is always positive, that 
is, Ea + E B > 0. This means that Ea > \Eb\, and thus only a part of the 
absorbed positive energy Ea can be retrieved from a black hole by this QET 
process. 

When we consider higher- dimensional spherical black holes, essentially 
the same results are obtained. The POVM measurement and operation for 
energy extraction are performed non-locally over a spherical-shell region sur- 
rounding an event horizon with a fixed radius. 




(31) 



/ takes a negative value, that is, 




(32) 
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4 Summary and Discussion 



In this paper, a new quantum mechanical method to extract energy from 
black holes with contracting horizons is proposed. This quantum mechanical 
method is based on a gedanken experiment on quantum energy teleportation 
(QET), which has been recently proposed in quantum information theory. 
Near-horizon regions of large-mass black holes, a POVM measurement de- 
fined by Eqs. (jSJ) and (jSJ) is performed for quantum fluctuation of N massless 
scalar fields in the Hartle-Hawking state, which is effectively described by the 
Minkowski vacuum state |0). Then, we obtain iV-bit information about the 
fluctuation denoted by D^. During the measurement, wave packets of the 
fields with positive energy given by Eq. (|T6|) form and then fall into the black 
hole. The amount of excitation energy is independent of the measurement 
results. After absorption of the wave packets and increase in the black hole 
mass, a D^- dependent local operation of the N fields given by Eq. ( |23p is 
performed outside the horizon. Then, accompanying the extraction of posi- 
tive energy from the fluctuation by the operation, wave packets of the fields 
with negative energy given by Eq. fl32|) form and then fall into the black 
hole. Because both the mass and the entropy of the black hole decrease 
due to absorption of the negative flux, this QET process is analogous to the 
Hawking process generated by pair creation of excitations with positive and 
negative energy, although the energy extracted by QET does not show the 
thermal properties. 

In general, black holes swallow various kinds of information from falling 
matter, which increases the area of the horizon. Because falling matter with 
the same mass, charge, and angular momentum generates the same geome- 
try in classical theory, classical black holes forget the information after the 
shifted horizon stabilizes. Meanwhile, many physicists believe, respecting 
the unitarity of theory, that quantum black holes do not forget the informa- 
tion and that the detailed memory is stored in some quantum mechanical 
ways. However, though many efforts to understand the quantum memory 
storage mechanism have been made so far, the complete resolution remains 
elusive. Even in the case of large-mass static black holes, the question as to 
where the information is stored in the black-hole spacetimes has not been 
completely resolved. The information might be stored inside the horizons, 
on the horizons, or even outside the horizons. This is the so-called black 
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hole entropy problem. From the results of this paper, we can conclude that 
some memories of absorbed quantum matter remains outside the horizon for 
a while, even after the shifted horizons are completely settled. The reason 
is following. If the unitary operation in Eq. ( 123]) corresponding to a dif- 
ferent iV-bit sequence D^'^k is performed after the black-hole absorption of 
wave packets with information D^, the amount of energy extracted by QET 
always decreases. Any quantum operation [/& with a wrong bit number b 
does not generate negative energy; instead, it generates positive energy of 
the field. Hence, this execution of Ub needs work. The amount of energy is 
E' B = 3t7 2 /4£(> 0). Figure 5 shows a schematic diagram for the case of a 
wrong iV-bit sequence = (Oil • • • 1). Thus, if the wrong bit number per 
N exceeds 1/4, there is no energy gain from the black hole on average by 
the QET process for the N fields. In an extreme random choice of Z7& 

yields rf /4£ energy input to the field on average. The black hole returns the 
maximum amount of energy by QET if the unitary operation obeying the 
correct D k is performed on the N fields. Therefore, the black hole remembers 
which Dk state was absorbed. This implies that the black holes absorbing 
quantum states with different D\~ are physically different even though their 
black-hole geometries are the same outside the horizons. In this case, the in- 
formation of Dk is imprinted in the quantum fluctuation of the fields outside 
the horizon. In fact, in the region out of the horizon, the one-point function 
of Il + depend on the measurement result D k as seen in Eq. ( 12"2"]) . Though this 
observation is of much interest, these memories may be, in a precise sense, 
not black hole hairs contributing to the black-hole entropy. This is because 
the memories are not stored eternally, but may decay in time T as suggested 
by Eq. ( 13 ip . This loss of memories over time suggests that the black-hole 
QET can be used as a probe for the thermal relaxation process of the black 
hole to a true quantum thermal equilibrium state. 

Before closing this discussion, I would like to add a comment about the 
Unruh effect in the Rindler spacetime [16] [17]. If we consider a detector 
coupled with quantum fluctuation of fields and uniformly accelerated in the 
vacuum, we observe thermal excitation in the detector as the Unruh effect. 
This effect is artificially caused by measurements and not a spontaneous ef- 
fect. The artificially induced effect is similar to the QET method described in 
this paper. However, the QET method is completely different from the Un- 
ruh effect. The detector of our POVM measurement is coupled with quantum 
fluctuation for a very short duration, though detectors used for observation 
of the Unruh effect must be coupled with fluctuation over a very long period 
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of time and should be simultaneously accelerated uniformly. The detector for 
QET is not accelerated at all. Moreover, the measurement results read out of 
the detector are essential for QET. We give a quantum feedback to the field 
based on the measurement result. Because of the dependence on the mea- 
surement result, the local operation given in Eq. ( 123]) can generate negative 
energy on average of the field with extraction of positive energy outside the 
horizon. However, the measurement result in the Unruh effect is not reused to 
give any feedback to the field. This is a quite different point from our model. 
In the reference [17] , Unruh and Wald also discussed a correlation of response 
of a particle detector with energy generated in a causally disconnected region 
from the detector of the Rindler spacetime. This causally disconnected re- 
gion corresponds to inside an eternal black hole when the detector is located 
outside the horizon. Interestingly, positive energy density is induced in the 
causally disconnected region when the particle detector outside the horizon 
is observed in the excited state. When the detector is observed in the ground 
state, negative energy is induced in the region such that the average energy 
density in the causally disconnected region remains zero which is equal to 
the vacuum-state value. However, it should be stressed that the energy in- 
duced by the detector is located inside the horizon and cannot be used for a 
person who stays outside the horizon. Moreover, the average energy (mass) 
of the black hole does not change at all after the out side-horizon measure- 
ment. This feature is completely different from that of our QET model. In 
the QET model, the person outside the horizon can effectively get positive 
energy from the black hole via quantum fluctuation of the field outside the 
horizon. The negative energy given by Eq. ( 132]) is generated outside the 
horizon by the measurement and falls into the black hole so as to decrease 
the black hole energy. 
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Figure Captions 



Figure 1: Spacetime diagram of the QET process in the near- horizon 
region. The future horizon stays at x~ = 0. At t — 0, quantum fluctuation 
of / in the spatial region [:ri,£2] is indirectly measured by a spin probe. 
A wave packet with positive energy Ea is generated by the measurement 
device; thereafter, it falls into the future horizon. At t — T, the unitary 
operation Ub, which is dependent on the measurement result b, is performed; 
this operation generates a wave packet with negative energy — I-EbI, which 
falls into the horizon by extracting positive energy from the field. 

Figure 2: Spacetime diagram of the horizon shift in the X coordinates. 
The initial horizon Hi stays at X~ = 0. By absorbing the positive energy 
flux generated by the measurement, the horizon moves outwards. Hp denotes 
the shifted horizon. After the absorption of negative energy flux generated 
by Ub, the horizon recedes. The shifted horizon is denoted by Hq ET . The 
positions of the horizons are fixed soon after the wave packets are swallowed. 
The argument does not change at all if the falling energy flux interacts with 
matter within the black hole. In the diagram, right-moving matter, which 
is represented by dashed lines, collides with the falling positive-energy flux 
inside the horizon. 

Figure 3: Schematic diagram of the measurement process of QET. The 
N cylinders in the figure represent the measurement devices used to detect 
the quantum fluctuation of the N fields. In this figure, the iV-bit output of 
the devices is D k = (101 •••0). Each detector creates a wave packet with 
positive energy Ea that is independent of b. The wave packets are absorbed 
by the horizon, which is represented by the shaded region to the left. 

Figure 4: Schematic diagram of a unitary operation process that is de- 
pendent on the measurement result. The N cubes in the figure represent 
the operation devices of the unitary operation for the N fields. After the 
operation, a wave packet with negative energy — \Eb\ is emitted from each 
device and is absorbed by the horizon. Before this, the horizon has already 
absorbed the wave packets corresponding to D k = (101 • • - 0) that are gen- 
erated by the measurement. Local energy conservation assures extraction of 
positive energy +|-Eb| from each field by the device. 
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Figure 5: Schematic diagram of a wrong operation on the quantum fields. 
In the figure, a unitary operation corresponding to an iV-bit data (Oil • • • 1) 
that is different from the correct data (101 • • - 0) is performed. Devices with 
wrong bit numbers do not generate negative energy flux, but generate positive 
energy flux with + |-E , #| = 3^ 2 /4^. When the wrong bit number per N 
exceeds 1/4, there is no energy gain from the black holes on average by the 
QET process. 
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